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The principal object of the calculations contained in the following paper, is to 
investigate the figure which a fluid, consisting of portions varying in density accord- 
ing to any given law, would assume, when every particle is acted upon by the attrac- 
tion of every other and by a centrifugal force arising from rotatory motion. To 
what extent this may have been the original condition of the earth, is a doubtful 
question ; and although observation does not fully warrant this supposition of the 
regular arrangement of parts, it has necessarily been made the foundation of most of 
the mathematical calculations connected with the investigation. Before proceeding 
to this problem, it is necessary to calculate the attraction of a body of any given 
figure, and consisting of strata, varying in their densities according to any given law; 
and it is in this problem that the principal difficulty lies. The elegant method of so- 
lution discovered by Laplace is well known ; and I have followed his steps as far as 
the point where the equation, known by his name, first appears. In order to illus- 
trate the nature of the deviation which I have there made, it will be necessary to 
mention some of the principal steps of the two methods. 

By means of a theorem, which Laplace laid down as true of all spheroids that 
differ but little from spheres, and the properties of the integral of the equation re- 
ferred to, he was enabled to substitute the easy rules of differentiation for the more 
complicated inverse processes, and thus to compute the attraction of that class of 
figures. It has, however, been since discovered by Mr. Ivory, that this theorem is 
true only of spheroids of a particular kind ; and, consequently, to this kind the solu- 
tion of the problem is restricted. This defect, and the indirectness of his analysis, 
led other mathematicians to consider the question; and, in 1811, Mr. Ivory pub- 
lished his method, which has the great advantage of being more direct, though 
equally limited. 

The method given in the following paper does not appear to be confined in its 
operation to any particular class of spheroids ; since the coefficients of the series, into 
which the required function is developed, are determined absolutely, without any re- 
ference to the form of the spheroid to which they are about to be applied. The prin- 
cipal change consists in the different manner of treating this partial differential equa- 
tion. Laplace and the subsequent writers on this equation, both as applied to the 

l2 
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calculation of attractions and the mathematical theory of electricity, suppose the co- 
efficients of every term of the series to be expanded into another series of the sines 
and cosines of multiple arcs ; and they avail themselves of the property which these 
terms possess of vanishing, in certain cases, when integrated between certain limits. 
The success of this plan, however, depends upon the restricting hypothesis above re- 
ferred to, that the radius vector of the surface of the body is capable of expansion in 
a series of terms, each of which satisfies Laplace's equation. The following method 
shows that the coefficient of the general term of the first series is independent of one 
of the variables, and thus dispenses with the second series of expansions. This re- 
sult I have arrived at, by first obtaining the integral of Laplace's equation in its 
most generalform, and deducing the arbitrary functions introduced therein, from 
considerations which enter previous to the equation of the surface of the attracting 
body. These coefficients being known, it is evident that the attraction of any homo- 
geneous body on a point within or without it may be immediately found when the 
equation of its surface is given, since it then depends only on a series of explicit and 
definite integrations of known functions, which can always be effected, at least ap- 
proximately. From this, the attraction of a heterogeneous body, similarly circum- 
stanced, may be found by the usual method of dividing it into concentric layers, and 
summing the several attractions of these, deduced as above. 

By Substituting the attraction so obtained, in the equation of equilibrium of a fluid 
body, Clairaut's theorem is immediately deduced ; and, from a peculiarity in the 
functions representing the attraction, it will be seen, that the same principles with 
longer processes may be carried on indefinitely, without the necessity of actually 
determining the precise form of those functions. 

The restricted species of spheroid above referred to, comprises all surfaces of revo- 
lution ; so that it is sufficiently extensive for most practical purposes ; but the inte- 
gration of Laplace's equation renders the analysis more direct, and the theory more 
complete. 

On the General Problem of Attractions. 

1 . Let f represent the density of a body at the point {x, y, z) ; and let /, g, h be 
the coordinates of a particle attracted by the body, parallel respectively to the axes 
x,y, z; then, if the power of attraction be inversely as the square of the distance, the 
resolved part of the attraction of the body, parallel to 

§{/— x)dxdydz 



V is err q{g-y)dxdydz 

y JJJ {(_/■_ ^)2 + (^ _ yY + (A _ ^)^}4' 

z is err qQi-z)dxdydz 
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the limits of integration being determined by the equation to the surface of the 
body. 

2. Let V represent the sum of the products of each particle by the reciprocal of its 
distance from the attracted point. 

Then V -JXf^ ^_^^^ +^(^l'^y)f + (^ _ ^)2)i ' and. by differentiating V, we obtain the 

d^Y d?Y (PY 
well-known property -jp + -j^ + -^ = 0, or — 4 t g', according as the attracted 

particle is not or is within the attracting mass ; f' being in the latter case the density 
of the attracted particle*. By transforming these equations to polar coordinates, we 
obtain 

^V , 2 dY , 1 d^Y , „ .^'^V , 1 ^V ^ ■ , 

and 

V — /*'■ 7"'^ Z*^ "- gr'^dt^ sine' dO' dip' _ 

"T/o^o »/o {r^+r'^ — 2rr' {coadcos 6'+ siad sin 6' cos{ip — 4>'))}^' 

where r^ ~p + g-^ + W,(^o% 6 — ,, ,.a g , ^gj , tan (p = y; and similar expressions 

are true of r' ^ and ?>' in terms of x, y, z. 
Put cos ^ = ft, and cos ^ = (ji/, and they become 

d^(rY) , d /,, o.dY\ , 1 <PV ^ . / 2* /, ^ 

/^r p\ pi.v qr^'^dr' dji.' d<^ 

^ —Jq J-\M {^s + ^2 - 2rr' (jn-ft' + v'(l - i>?) -/(l - /*'^) cos (p - ip'))}* ' ^^'^ 

3. Expansion by the binomial theorem shows that 

{^2 _j_ ^2 _ 2 r/ (|«-/a' + -/(l — (M.2) V {I — |«-'2) cos (?) - ?)'))}* 
may be expressed either in powers of r or of/; thus 

Po^ + Pi;J + ..- + P„;7H:i + --vOrPo- + Pi^ + P2;5 + ..-|-P„^:^ + ..., 

where P„ is a symmetrical function of ]«/, \^{l — i^"^) cos p, '/{I — |«-2) sin (p on the one 
part, and f/j', \/(l - /*'2) cos ^', -/(l — ,!a'2) sin ®' on the other. 

By substituting the first expansion in (2.), and the value of V so obtained in (1.}, 
we have a series of equations 

= 0, or — 4 5rf', 
except when « = 2 ; and in all cases 

which is the equation of Laplace's coefficients. 

* Vide Pratt. Mec. Phil., § 168. Laplace, M^c. C61. liv. iii. t Vide Pratt. Mec. Phil., § 169. 
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4. This equation was not integrated ; but by a skilful use of its properties, the 
problem of attractions was greatly simplified by Laplace. He laid down a theorem, 
respecting the surfaces of all spheroids of small deviation, that their radii vectores 
might be developed into series, every term of which would satisfy the above equa- 
tion ; and he also gave a method of expansion. By means of this theorem, the pro- 
blem could be solved for spheroidal bodies which differ but little from spheres ; but 
its generality has been greatly restricted by the researches of subsequent writers*, by 
whom it has been shown that it is true only for bodies whose radii are expressible in 
rational and integral functions of (t},\^{\ — fo'^) cos ^', Vi\ — (m'^) sin p'. Among these 
are the ellipsoid and elliptical spheroid, and a large class of other spheroids. In these 
papers I have adopted a different proceeding ; I integrate the equation itself generally, 
and determine the arbitrary functions contained in the integrals by the circumstances 
of the problem itself. In consequence of the peculiar form which P„ then takes, V 
may be found by effecting the operations indicated, which are only explicit integra- 
tions. 

5. I shall now proceed to integrate this equation. 

Consider [a and p as functions of two new variables X and Y, to be determined 
from the equations, 

dX = ^(dp + kdf.) 



(If 

Y 

lip 



dY = ^{dp + l^di.), 



where k and k' are the roots of the equation (1 — u.^) ^2 _j_ ^ _. q These roots 

1 — /A 



are -= — —J^, whence we obtain 

1 ~ !"■ 



X = ^ + iN/-lHH7,andY = ?>-ix/-llogf^. . . (4.) 

d(>,^ — \d[t.J dX^ "^ '^ di>. di>. dXdY "<" \diJ.J dY^ "'" dp.^ dX ^ dp,^ dY' 

_ _i ^P„. ■ 2 ^P. 1 c^^P 2^>/1T /rfP„ dP\ 

- ~ (1 - p.J dX^ "^ {I- y,^f dXdY (1 - p.^f dY^ "^ (1 - t^'^f \dX ~~ Iy) ' 

dV„ dXdP„ ,dYdP„ _ i/^n[ (dP„ dF\ 
dix, — d[^ dX^ dfu dY ~ l-i>? \dX rfY/' 

rf^p„, (dxy^v dxdY d'p„ (dYyd^p d^xdP d^Ydp„ 

df —\d<p/ dX^ '^ '' df d'p dXdY'T' \dfj dY^ '^ df dX "•" df dY' 

— dX^'^'^dXdY'^ dY^ 
Substituting these in (3.), we obtain 

4 d^F 

* See two articles by Mr. Ivory in the Philosophical Traftsactions, 1812. 



AND THE FIGURE OP THE EARTH. 79 

From (4.), by subtraction, 

whence 

_ s-(X-Y)V^l_l 

and 

1_ 2_ 4.-(X-Y)V— 4 1 

2 

consequently 

J-V J AT "1" 

4 9 ' 



S$^+"7^^^ = «- Letn (»+!)= a. 



Let 



or 

Let 
Then 

and 

whence 



2 



dY = '^^ <^ben^ + J_Y = 0> and P„ = - ^ - cos^ -^- ^ 

4C0S _^— 

<iP„ tZ^p 4 oX-Y dv 4 X-Y . X-Y 

^ = dY = - dXdY -^^OS -2 dX-^ cos — ^sin -2—' 

d^v dv^ X — Y at) 

dXdY + dX ^^'^ ~~2 T , ,X-Y — "• 



4 cos^ : 



do , . X-Y ^ 
^ + t;tan-2— = #. 

d^v dv X — Y y df 

dXdY + dX^^ ~T~ + ";;^ ^X=:Y = dX' 

2 CDS'* 

2 

dt . a-2 w ^ 4 „X-Ydi. 

_ + 4 _,x~Y = 0; orw = _ ^-^cos2-^^. 



cos^' 



2 



or 



Let 



d7) 

dY — 


4 
-«-2« 


t= - 


a -2 COS 


d2^ 


.4. 1^.9t 



2X-Y_^^ 4 X-Y . X-Y d^ 

^s 2 dXdY~a-2C0S-^— sm— 2— ^; 

X-Y d'H 8 X-Y . X-Y dt 

— 2- dXdY - o-^Ts COS -^- sm -2- Jx' 



X-Y , , a-2 _ 

dXd Y "T dX ^ """ 2 + * "; aX— Y - >^- 

4 cos^ 



d; , o/* X-Y 

_ + 2nan-2- = 9. 
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and by repeating a similar process, we obtain 

7 

4cos^ 



d'q . dq X-Y a-6 

dXdY + dX^^^^^2^ + 9 , ,X-Y-"- 



2 

By observing the assumptions here made, and the results obtained, we find that in 

/dV \ X Y 

the first assumption y-^ = v), the coeflBcient of P^ tan — ^ — ^^ ^» ^^ *^® second, 

V Y 

that of V tan — ^ — = 1 ; and so on, in the order of the natural numbers ; and in the 

results, the numerical coefficients are 1 , x > 2, —7 — ; 3, — j— . . . generally^, 7 -. 

I shall prove this in the general case, by showing that if it is true of one value of n 
(as we see it is), it is true of the next value, and so on. Let the {n — l)th substitu- 
tion give 

-gi--Lig-fa-lUan^^ + . ""^"~^^^"~'^ = 

4 cos — X — 



and let 



then, as before, 



d§ , , :. ^ X-Y 

^ -f (n - 1) f tan —j- = si 



d^g , . ,,^ X— Y do , n-\ 1 As. 



and, therefore, 

A + r 7x-Y = 0, and ^ = - ii^zrT^inw cos^ 



ds , a-{n-l)n ^ , ds 4 ,X-Y, 

— " , X-Y = Q^ ^nd g = -Jx a-in-l), ""°' 

4 cos^ ^ ' 



, X-Y . X-Y 

4 cos — ;: — sm 



dg _ d^s 4 , X-Y ds 2 2 



■cos 



dY — dXdY a -{71-1)71^^" 2 dX a-{n-\)n 



Consequently 



,X-Y ds 4 X-Y . X-Y 

cos^ — jr- Tv ' 7 rrr cos — -x — sin — 7i — 



or 



dXdYa~{n-\)n^^'' 2 dXa-(n-l)n 

ds 4(«-]) X-Y . X-Y _ 

^ dXa-[n-l)n^^^ 2 ^'^ 2 ~ *' 



d^s 4 oX-Y . ds 4« X-Y . X-Y 



cos'--o— + jvTir7r-:TwC0S--H--sm— 5— + *= O; 



dXdYa-{n-l)n^'''' 2 T^rfXa-(w-l)M 
that is, 

d^s , (fs ^ X-Y , a-n(n-l) „X-Y ^ 
jX^ + ^«tan-2-+*— -^ cos2--2- = 0, 

and, therefore, the law of coefficients, as above stated, is correct. 



Flul.^ r^MximCCCXLL Flute IT,, p. Si. 



5^ 



^ 









^ 
^ 
S 
^ 









'^ 



■^ 






I 



SI. 

i 







1 

■13 

I' 



y/.Boyt^rnmb ad/,na£dA 
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Restoring the value of a, we get 

d^s ds X-Y , s.n 1 

COS'' 

^g X Y 

which is integrable ; and its integral is ^ 4- w . * tan — § — = some arbitrary func- 
tion of Y, as X Y. Integrate again, then 

- fnt 

where -^ is arbitrary. Effecting these integrations, and reducing, 

« = cos — 2— (^y cos -— 2— .%Y(?Y4.^|^Xj. 

To return to P^, we have the following systems of equations : 
V=fvdY, 

Wrrrg •/ 2 \ftiJ 2 rfY^ = COS —^—JtCOS^—j-dY, 



itiai'^~^ 






#= COS 



_4Y-X/' .Y-X 



2 



- / 9 cos* — 2 — . rf Y, 



whence 



^^^_2(„-i)Y-Xy>^^^^2(»-i)Y-X_^Y, 
cos-2»^(/cos'^»^.^YiY + ^x); 



P„ = .....ycos-2^^^cos-2^5^(y^os2»^^5iY</Y + 4x)<;YiY...(« times.) 

Now cos ^ = cos (y^n log \/4^) = I (v^-f^ + \/ttS)' 

Y X 1 

and cos^ „ = . _ a ; and the complete integral will be expressed, by substituting 

for X and Y in terms of jO. and (p. 

6. But an important point yet remains to be determined. The original equation, 
being a partial differential equation of the second order, can only involve in its in- 
tegral two arbitrary functions. But here, after % Y and -^p X have come in by two inte- 
grations, we have n integrations to perform with respect to Y. It would seem, there- 
fore, that no constant or arbitrary fimction of X must be added in these integrations. 

MDCCCXLI. M 
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Such is not the case. At each integration a function of X must be added, and these 
functions determined by reference to the original differential equation*. 
7. Returning to the value of P^, we have 

Po = %i {p -j^~^ % r^) +^{p + j^-^i log 'rH)- 

Now in the calculation of attractions, where Pq is the coefficient of r^ in the ex- 
pansion of 

{r2 ^fJ'i ~rr' {y.^' -\^ ^ (\ — ^i) ^ (1 _ ^'2) cos {p - ?>')) } ~ *' 
we know that it is 1 ; consequently 

X, (* - |n/^ log [^) + -^ (p + |V^ log 14^) = 1, 
and expanding by Taylor's theorem, we get 

+&v-=iaiogH^-:y+a.. 

By equating the coefficients of the same powers of g- log , _^ , we have %i <p + 1^ ?> 

= 1 , and 4' ?> — %'i <? = 0, or •*!/ (p — ^(^j (p = constant. 
Therefore •»!/ (p and xi *P are absolute constants, and their sum is 1 ; whence it fol- 



^ = 1. 



lows that ;^ (p = 0. Let -^ <t> = », then 



- = K • • • -y ^o^" 



2Y-X 



- / cos" 



.2Y-X 



^ 



cos 



-aY-X 



<ZYt?Y . . . (» times)). 



2 J ^"" 2 c/. -2 ^^" 2 

Effecting these integrations, we find that P„ consists of a series of powers of 

* The common differential equation (1 — ju«) ^^ — 2u^^ + « (» + 1) P« = will, illustrate this point. 

dfi^ a fi 

Let Pm = -j -, and after substitution, differentiate n times ; then (1 — /*^) — - — 2 (» + 1) u — - = 0, whence 

a fi—n dfx^ d IX 

z = / ^ , , + m. It is clear that no more arbitrary constants than k and m can be introduced ; and 

yet if the integrals were left indefinite, we might obtain an integral of an expression which should differ from 
the integral of the same expression obtained by a slightly different process, by a constant. By another inte- 
gration this would cease to be a constant, and we should obtain thus different values for ?«. The fact is, that 
constants must be added at each integration, and recourse had to the origmal equation, to determine them in 
terms of m, k, and ^i. 
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Y — X 

i — g — , whose coefficients may, for anything we yet know, be functions of X. 

_i ^!-_. ^__.L coefficients will satisfy the original 



tan — ^ — , wuose uoemuienis may, lor anyining we yei Know 

The following process shows that constants as coefficients will satisty tfte original 
equation, and determines them. The integration itself gives the coefficient of the 
highest power. 
Let then* 



* 

tan" ^^ + d tan» -^ ^^-^ + c" tan« ~^ ^-^ 1 
I +... + c(»-^)tan2^^^4-c(»-'>tan^^ + c«J* 






r 2»-2 



dV. 



ilY 



^ = c 



tan'^-^^^^ + ^'^'tan' 



ljt^^n-2^L:zA 



-^ ... ,-3]L^+....+ |e(»-)tan^[(l + tan^^), 



^c"tan» 



^ 2 



+ -C 



1 Jr^ - 1) 
2 



'^l::^tan»+^'^? + ^ 



= c ■ 



+ 



+ (^'c"' + ™-^c')tan» 
(l^,.+»-:z_%.)tan«-3]L= 

"^ - ~- ' («-^))tan3 



^-i^<^'tan»^+e-i-%" + ^)tan»- 



Y-X 



2 



Y-X 



i.c 



+ 

Y-X 



Y-X 



Y-X 



..+ (|c(»-3)+|c(»-^))tan*^ 



+ a^ 



^» - »)) tan2 ^5^ +1 c(» - 2) tan ^^ + i c(» - ^) 



,(«-i)+|cO -;.«» 2 






c 

Y^^ 



dtSiXl^-^^-^ 



'(iLp)^ tan» "Lll + 'ii^) ,- tan' 
[re — 1] 2 ' 4 

+ -T- V-^ ^ + F=ri]j tan» 



- 2 /« - 3 ^,„ 



2 
Y-X 



(^c^" + ^c')tan»-3^ 



^ ( lc(« - 3) + |c(» - 5)j tan3 ^^ 
# n«i = 1 o 9 J. « 



n o o ^ 



[«] = 1 . 2 . 3 . 4 re. 

M 2 
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M+ 1 n 

n — ;; — P 



« aY-X 



COS' 



Y-X 



(^+l)^%an'' ^ + '^-^ c' tan»- ^^' 






...+ »_(!i+i) e(-i) taa^ + "^) cW 



Whence we obtain c* = 0, c"' = 0, c' = 0, &c. 
Also 



/w(w+l) (w -l)(w-2) \ „ : j^--3^ / ra(ra+l) (w-3)(w-4) \ j^ _ (w-3)(w-2) „ 
V4 4 /^— rw-2VV4~ 4 ) ^ — 4 ^' 



( , 

Consequently 



[m-2] 

w (w + 1) _ (w — 5) (w — 6) 
4 



).- 



(w - 5) (w — 4) 



c'^ &e.&c. 



^^"~\an'^^:;^ + 



Now 



2»- 



-1 



+ 



2 ' 2[w-2](2m-1) 



tan' 



»-2 



Y-X 



2.4. [w - 4] (2 re — 1) (2w - 3) 



tan 



K — 4 



2 
Y-X 



2» - 1 tan' 



«-6 Y-X 



Y-X 



tan 



Y-X 



+ 2.4.6.[ra-6].(2w-l)(2w-3).(2re-5)' + ' 

= -iN/^^log5-^= tan-1 (-^^^-f),... 

whence, finally, 

^ W 2. [re-2].(2w-l) "^2.4. [re-4](2re-l) (2 w-3) "•"■■■/' 

a remarkable result, showing that in this instance P^ is independent of p. 

P^ being free from ip, is a perfectly symmetrical function of jo, and |M,' ; and ,«,' is a 
constant with respect to /«- ; therefore 

^»— *%V [„] ^2[re-2](2re-l)^2.4.[w-4] (2m-1) (2m-3)^"V -^ 

V [«] ^ 2.[k-2](2«-1)^*"'/- 

To determine K^ for any particular value of n, we refer to the expression from 
which the two series were deduced ; namely, 

|r2 4- r'2 — 2 rr' (jotju,' + ^/(l — ^i?) a/(1 - j«,'2) cos (?) — ip')) } ~ *' 

When j«- and /*' are each 1, then P^ = 1, which gives an equation to fiind K^. 
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8. Returning to the theory of attractions, we have, when the particle is internal, 

= 2T/''/'^f{Por' + Pir + P2? + .. + P„;?^ + ..)^r'^f*', 

for that portion of the body which is comprised between a sphere of radius r, and the 
surface of the body, supposed to be a surface of revolution round the axis of 2. ; for, 
in that case R, the value of r' at the surface, is independent of <p'. 

9. Suppose, for example, we wish to find the attraction of a homogeneous spheroid, 
on a point within it. In this case § is constant, and 

a being the semi-major, and c the semi-minor axis. 

First, all the even terms vanish ; for the general even term is 

/I — U.'* (2\ 2»-l 

Now Pgn . 1 consists of odd powers of fJtJ ; and y — ^ 1- ^j 2 can be expanded 

in even powers of (jtJ ; therefore the integral of the product (which is an odd function), 

taken from [jt>' — — I to (J = 1, is 0. Also / _ ^ Pg^ , j dfJ = *. All the odd terms 
above the third vanish ; for the (2 w + 1 )th term is 

+ iT^y^ 1 P2«''^ ^i^'' f'*'" ^* '"^y ^^ shown that/"^ ^J'^'' P, P.' dfjb' dp' = 0, if < and / 
be different integers. Now when n is greater than I, \-~^— + -^) is a rational 



and entire function of (jij', and, therefore, capable of being expressed in a series of 
Laplace's coefficients -f-, the highest of which will be of the (2 n — 2)th order ; and 
therefore no term of this expansion can be of the same order as Pg „ ; and the inte- 
gral of the product of any two of different orders, between these limits, vanishes. So 
the second member of this vanishes. 

* See Pratt. Mec. Phil., § 180. 

t See Pratt. Mec. Phil., § 176. Poisson, Thferie Math, de la Chal., chap. viii. Laplace, M6c. C^l. hv. iii. 
chap. ii. 
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The first term is 

— 2'!r§ia^ ^^ — ry, {e being the eccentricity). 

The third term is 

t«„-l -. . _ -), 

= -3t^(|-^2) ^2 (i-.-^ + ^^iL^sm-^e), for K^ = 9. 

4: Tt P y*^ 

The value of V for the sphere of radius r, calculated by the usual method, is — j — ; 
consequently, for the whole ellipsoid, the value of V is 

-^ +2r§[a^-^^ — ^sm-ie-r2j-3Tf(^y-fi2j^2^____^^___^sin-iej. 

Differentiate to/, by means of the equations r^ =f^ + g^^ + ^^ and ,0, = — and we 
have 



= - '^ ? f (t - ^0 ^' (i (-^r^ tan ^ y ^::^, 9 



— -w. = attraction in .r = + 2 "^gfil — — H -~j — sin~ ^ e) ; 

ion in 3/ = + 4 cr ^ A (+ -^ -^ sin- i e) ; 



df 

dY . . 

SO — -7— = attraction m 

dg 

so — -j^ = attraction 
Also 



which are the 
common ex- 
pressions 
otherwise 
found -^. 



10. By a similar process, I have deduced the attraction of an oblate spheroid, on 
a point within it ; the density varying inversely as the distance from the centre. The 
corresponding expressions are 



dV 



J 3 ^§/ 



(//— "f'^'^^ r ~ 4 a(l-e2)| e^ ~ 3 ~ V 2e^ "'"3 e"/'^^l-ej' 



~dff — '^'^^^r~4a{l-e'')le^~S~\ 2^ +3 e /'*^&l-e/' 
* See Pratt. Mec. Phil., § 172. t Ibid. § 158. 
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rfV , 3 TTg (3/^^-1) f 1 1 /(1-eT , 2 1 - e^\ , l+e1 

-^=+2t| + -4 -^(1 _ ,.) I ^2- - 3- - (^-^^r- + y —^) log fzr; J- • 

11. When the particle attracted is external, the series in (3.) does not give a finite 
expression. Instead of taking it separately, I make it a case of a general theorem 
which follows. 

12. To find the attraction of a spheroid on a point within it, the density being any 
function of the distance from the centre, and the eccentricity being small. 

Let f ?i r' represent the law of density ; then the value of V, for the portion com- 
prised between the surface and a sphere of radius r, is 

Integrating by parts we have 

/^r^/-^rfr' = ^ + («^l)^tf^-(«-l)n^i^ + («-l)n(«+l)-5^ 
Therefore the (n + l)th term of V is 

NowR = 1^ + (^^-i) ^'2 j ~* = a (1 + e>'2)"*= fl (1 - -|^^'2), rejecting e* 

and higher powers of e; and <p^^^ R = ip^^^ a — ^(^_X) « • "^^-^ to the same degree of 
accuracy. The last member of the expression for V need only be calculated when 
w = ; for all the rest of the terms (^involving /_ j P„ fi^i"-' where n => o) vanish. 

The first member need only be calculated when w = 0, and when ra = 2 ; for when 
n is odd, it vanishes as before ; and also when n is even and greater than 2 : for the 
functions of R involve no higher powers of \j} than the square ; and consequently 
they vanish, when multiplied by P4, Pg, . . &c., and integrated with respect to ^, from 
— 1 to + 1*. 

When « = 0, the term is 

2 -^ if_ 1 (R ?, R - <p, R) d(iJ - 2 ^f 7^1 {r<p^r - (p^,r) d(^', 

* See Art. 9. 
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When n = 2, the term is 




To calculate V for the sphere whose radius is r. 
Let the sphere be referred to polar coordinates, the centre 
„ being the pole. 

LetOC = r,OP = ri,andPOC=0; thenPC= ^/(r^-\.r^2 
— 2rriCOs6). 
Mass of the element at V = §ri^ pridviSinOdd d a, and 

_ 2v§f' '^' ((r + ri) - {r-rS) dr^ = ^X' r^^pr^dr^ 

K being the value of r^ ^^ r^ — 2 rj p,^ rj + 2 <p,„ r^, when r^ = 0. 
The whole value of V then is 

4 T^ ja^^a - «^„« - -|-e2a2<|)a - [(^.r - |-<|),,r + -) - ^ r^ (/*^ - |) e^ <|) a j . 
And the attraction toward the centre 

where 

•v|/ r = 2 </)jj, r — 2 r ^;, r + ^^ ^/ J' — K. =yo r^^rdr. 

13. To find the attraction of the same spheroid on a particle without it. 
The series (Art, 3.) is 

Now 

^f<l> r' .r"' + Kdr' = r'» + i <^^ r' - {n + 1) r'" <|)„r-' + « («. + 1) r-'» " ' ^,„ / - &c. (A) ; 

and the general (reth) term of V is 

^/_\P„_i (r'» + ' <P/ - (« + 1) r'"0,r' + « (» + 1) r'«-i <|,„,r' - &c.) d fyJ 
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where C is the value of (A) when r' = 0. As before, the last term need be calculated 
only when w = 1 ; and the first when « = 1, and « = 3. 
When « = 1, it is 

When » = 3, it is 
^{l.^-^)'^f_\df.' {i.'^-^) (R*ct>,R-AR^4>„R+ i2R^4>„R-24R<Pi,R 

The whole value of V is 

^l^; fl - "I" a3 <|, a _ - (fi2 _ _.j _^ a . e2 j. 

And 

14. Instead of the eccentricity e = y^ ( ^ — ^) i* will be more convenient to em- 

c ^ 

ploy the ellipticity s = 1 — -*. These give s = y' ^'^'^ *^^ values of V become for 
an internal point, 

Ax^>^a(^,a-<^,,a- (<^;, r - | (/)„,r + y) - -| a^ <^ a - ^ (^2 _ _j g ^a j . 
for an external point, 

15. To find the attraction on the supposition that the body is composed of sphe- 
roidal layers, homogeneous in themselves, but differing from one another in density 
and ellipticity. 

First, on an internal point. 

Let r*, as before, be the radius vector of any layer ; d its equatorial radius ; § <j> a' 
its density, and g' its ellipticity, being some function of a' as % a'. Then 

a' = r' (1 + g' |M-'2) and (p a' = <(> r' + r' (}>' r' x r' jm.'^ = ^ r' + F / . (m,'^^ suppose. 
Consequently to the term before produced in (12.) by ?i r' we must add a term 
similarly produced by F r' . (jJ"^. Also, instead of taking, in the first instance, the por- 
tion comprised between the surface and a sphere of radius r, we must take the por- 

* See Puissant, vol. i. p. 259, where the word ellipticity is used in this sense. 
MDCCCXLI. N 
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tion between the surface and that spheroid on which the point lies, whose ellipticity is Sj. 
The (n + l)th term of V now becomes 

a f, a - f„a — J a2 ^ a - y (fi2 — y) r2 ^ a 
-(a?)^a-«?)„a - ^ a? p a — -^ (i^.'^ - -^y r^ p a) 



The first part of this gives 4irg 



a being semiaxis major of the stratum on which r lies. To determine the other part, 
it is necessary to compute it when « = and n = 2, which gives 

To the sum of these we must add the value of V for the inner spheroid ; and for this 
purpose we have to obtain V for an external point. 
To the expression in (13.) we must add 

to be calculated when n— I and n = 3. This is 

The whole value of V is 

•4/a - —a'^pa- y(ju,2 — ^p^P a + -j f^ ¥ a! .a'^ d a! 



47r I 



H 



>. 



r--^ 



+ \-~;^fyo;,c^*da! 



After writing a for a, and Sj for g, add this to the other value of V, and apply the equa- 
tions 

f^^d . a! da! =Jl''a"^ s' dp a' = a'^spa- aH^pa -f^P d d (a'2 s'),, 

-Jl'^d .d'^dd = \f^" dUdpd^aU^pa-^\f^\dd{dU), 
and similar equations for the other integrals ; and we shall obtain 

apfl-p,,a~ap,a-^p,,a-\-^-^-\^f^\dd{d'^i) ^-\ f^\dd{dN)\ 1 
-T{^^~^{f.^^^d^-^^f,^^dd{d^i) } J 



V = 45rf 
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16. To find the equation of equilibrium of a heterogeneous spheroidal mass of fluid, 
revolving about its axis, with an angular velocity <w. 

By the principles of hydrostatics the general equation is J -r = / Q^dx -\-Y dy 

+ Z dz), ^ being the density at the point {oc,y, %), p the pressure, and X, Y, Z the 

sums of the resolved parts of the forces ; which are — ^, — -j-, — ~r-, and the 

centrifugal force. Let the axis of z be that of rotation ; then the centrifugal force is 
ft)2 X along X, and a^ y along y. Let us express a in terras of the ratio of the centri- 
fugal force at the equator to the equatorial gravity. Call this ratio m, which is small 
in the case of the earth, being of the same order as i. Then 



2 45rgm4/a 

mass 4 TT g 4> a' cfi 



m, = rrrr: = j-tt -rr, or &>'■ = g 



Therefore 






and 

^owj -J- is a constant for a level surface. Hence for any stratum we have 

C _ V -I ^^^sr^^ (1 — i^V r ■ 

At the surface this is 

= f?M«-5V(("-¥)N« + jf$Ma(l-K), 
where 

Ma=y^ <^a'-^^^^^ia',andNa = y^ c^a'-^'da'. 

For r write a (1 — s fA^), then 

Equate the coefficients of iJt>^, then 

a V 2 J — 5a^ ^^'*- 

17. By differentiating and changing the sign of y -r, w^e obtain the amount of 

gravity which acts towards the centre ; which, to the order we are now considering, 
is the same as the whole force of gravity ; since the cosine of the angle of the verti- 
cal differs from unity only by terms of a higher order. 

n2 
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Consequently 

CMa 3/„ l\Na m ,, ,. rMal 

TMa ,, . ^ „, 3 / „ 1 \ N a m ,, „, Ma") 

= 4 Tf |3^ (1 + 2 3^2) - y (^2 - -^) ^ - y (1 - ^2) -^ j, 

where G is the equatorial gravity and / the latitude. 

18. From this it appears that up to terms of the 1st order, R = a (1 — s sin^ /) is the 
equation of the curve which generates the surface of equilibrium, where the value of 
B depends on m, or on the velocity of rotation : but as the coefficients of higher 
powers of sin / may be considerable, it will be useful to find the surface of equilibrium 
to a greater degree of exactness. For this purpose we must introduce the fourth 
power of sin /, whose coefficient will be of the second order. Let the equation of the 
strata be «' = r' (1 + «' [i>'^ + A' ju.*), or r' = a' (l - 8>2 + (g'2 _ A') (U,*), where s' and A' 
are functions of a' as % a' and 6 a'. Then 

<|)a' = ^r' + Fr'|!A2 + r2-^2^^j^^ Y + ^^'^0»"i^^ = ^^+F»'-it*^4-nr./M,*,suppose, 

The (« + l)th term of V for an internal point now becomes 

By writing for r and R their values in terms of a and a, it is easily found that 
/'^i, d r' equals./;- ^,a-+{^){~. ,'' + p« (." - A') - ^■' "-^^ .') 

+ ( n-s ) 2^ -g — the same functions of a and Sj. 

And that similar equations are true for the two remaining integrals. 
1st. Let » = 0, and we have 



'^sj-i di^'- 



y^\' ((> a' d a' + a^ <p a[- e fji,'^ -\- [J* (j B^- A)) + a^ <j)' a' 



V'"' 



2 

+ iju'^/^" a' . F a' </a' + /!<,'2 a2 F a (- s fjt,'^) + v^^ f^ d "Hd dd 

— same function of a and gj, 
which gives 

4fff |y^''a'?>a'<^a'+a2^a(-y + y(|-g2_A))+a3^'a^ + |-^°Fa'.a'.rfa' 
— a2 p a y -f —J a' . n a' (^ «' > — 4 ff f (same function of aj Ej and Aj). 
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2nd. Let n = 2, and we get 



which gives 

*'l {l^ -l) '''{l"'{-i+k(i~ ^)) + ''f'''k'' +^fo' ^'''•' - k^ " ■ • 



+ ' 



„r / —J- da' I — 4 TT §\jjfi — Y ) r^ (same function of a, Sj and Aj). 
3rd. Let w = 3, and we have 



or 

— (|a,4 — — 1^,2 _j_ j 4 TT^r* (same function of a, Sj and Aj). 

19. To the sum of these must be added V for the inner spheroid, for which we shall 
have to find the value of V generally, when the attracted point is without the body. 
The general term of the series for V is 



r 

and 



f^r"'(pr' dr' =X"^'" •?>«'• <^«' + «"+ ^ ?=«(- ef*^ + f** (^^ - A) + [^'^y'^) 

+ a- + 'p'a'^; 



and similar equations hold for the other functions. 
1st. Let n = 2, then we have 

10 



'■^y- id(^' \f^" a'^ ^ a! da' + a^(p a(- i 1/.'^ + ^J^ {2i^ - k)) + a^p' a- 
+ y V2 F a' 6?fl' + a3 Fa (- s^'2) ^fj a'^ n a' rfa' 1 , 



which is equal to 

^ I ^\'2(pa' rfa' + «3?5a (- y + y (2^2 ~ A))+ a*?5' a^ + yy%'2 Fa' . <^a' 



- |-a3Fa.8 + yy"'^'^n'*''^^'} 
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2nd. Let n = 4, and the term is 
I- (^^ ~ y) ^^/V^' {(^''-^)[X''a'^<pa'da'-^a^<pa(-Bi.'^-\-i^'^(3B^-A)) 



or 

/U.2 — 



- gg «5 F a . s + 35^" a'4 n a' </a' I . 
3rd. Let n = 6, and we get 

+ a7(f>a(~ gp'2 + («,'4 (4 s2 _ A)) + a^ <p' « !!|i! + J'^" a'<^ Y a' d a' - a^Faaf*'^ 

which is 

20. In these three terms write a for a, s^ for «, and A^ for A, and add them to the 
other value of V, and apply the equations 

^"«'2 na'da' = \fj' a'4 d {% ^'^ <p' «') -^ f^^ a'3 ^ a' <Z (?) d) ~ a^e- (d a ■\- a^ kip a 
— ^fj d^x(^'^-d<pa! —J'J <pd d{d^ M) = —a'^s^p'a + a^ A(p a — 2 a^s"^ <p a 

+ 2^" (p a' ^ (a'3 2'2) -f^^dd (a'3 A'), 

and similar equations for the other integrals, and the value of V for an internal point 
becomes 

_e'+|-(2£'2-A')l) 

^.-Jldd 



47rf 



J^ ?«' L- j^ -^dd^^J^ pd -L-- ^ 

= Air Aw a- M' a + ^;. M a + |- (fi^ ~ i-) {r'^ N' « -™ r^ N' a + -i N a) 
where the functions M, N, P, &c. are substituted for the corresponding integrals. 
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21. The equation of eqiiilibriiim will be formed exactly as in (16.), except that 
the expression for the velocity must be found more accurately. By (17'), the force 

of gravity at the equator is 4 tt ^ 5^2 ( 1 -j- s ^); and the centrifugal force is jy^ a ; 

therefore m = jj— r 51^, and «2 = 4 ^ ^ ^-^ {^rn + m b ~ -^j. 

The equation of equilibrium becomes 

-^ = C = V4-2T^r2(l ^i,2j-^^m[l+B--^), 
which at the surface becomes 



3r 



1 / IxNffl 1 / 6 2 , 3\ Pa , 1 . .. .Maf , 3m^\ 



= f|a+^P' + A.^^)-i(^^-|)(l+3sf.2)^ + |(^4_|^2 + __, 



3\Pa 
35/ a^ 



+ "2— {I — 2 B (M^) ^s \^m + m S — -y-j. 



Equate the coefficients of fjb^ and those of ^*, and we have two equations for deter- 
mining s and A, thus showing that equilibrium is possible. These are, 
Ma/ m 3 ,3,\ Na/1 s\ 2 Pa 

, Ma, . , , Na 3s , 1 Pa ^ 

and i^(A + »i2)--^"--5-+ 9-^=0; 

, Ma/ m , 6A , 3 , , Wis H o\ Na 

whence ^ f^s _ -^ +_-^ + -^ ^2 4. __ _ _. g2| = __, 

M« /a o 9 1 5ms\ 1 Pa 

22. The resultant attraction in the direction of r, is obtained by differentiating 
/— r, as found in (21.), and changing its sign. This produces 
^ fMa 3/ l\Na 5/, 6„,3\Pa ., „, / , 3m^\Ma1 

since the terms arising from differentiating the expression for V in (20.) with respect 

to a, cancel each other. 

For r write its value, and arrange the result according to powers of (*, and this 

formula becomes 

"Ma/, , 3»i® , „/^ , , 3m*\ , , .^ . , o \\ 

3-^2 \l —m-\- ms g- + (Jb^{2i + m -\-2ms s- > -\- [i* {2 A -\- s^ — mi)) 

Na/ 1. 2^3 A ^ 12\Pa/l o^^ , . 5\ 

Maj- 37» 27 . 9 2 3 

= 4 Tf j^^l 1 + s - ^ - j^ms -f 4-W2 __g2 + _, A 



47rf 
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by applying the two equations deduced in last article. The total attraction normal 
to the surface, or the force of gravity, is found from this by dividing by the cosine of 
the angle of the vertical, or by 1 — 2 s^ (|«,2 _ ^4y . whence 

Mar,, 3m 27 ,9, 2,3. 

Let K be the sine of the real latitude ; then as fi^ = K"^ — 4 e k^ — X*, we get 

Ma r 3m 2? , 9 o 2 2,3 

g- = 4 T^ 3"^ |1 + g - -2" - Yj m g + ^ m2 - y s2 4. _ A 

Let G represent the equatorial gravity ; then 
g- = G|l + X2(|m-s--yma + ys2 + yA)-X4(2e2H-3A-|-ms)| 

= G|l + sin2/(|m-g- J|»ie + yg2_|-A) + sin2Zcos2/(2s24.3A-|ms)|, 

which is an extension of Clairaut's theorem. 

23. In this process A indicates the amount of deviation of the required surface 
from the surface represented by r = a (I — s fjb^). If the equation had been taken 

a=:r< 1 + e|M<2 — (^ — e2 + Bj (jm,^ _ ^2 j _j_ £2 ^4 I g would have been the index of 

deviation from an elliptic spheroid. 

To apply Clairaut's theorem to this surface, we have 

3 e'^ 

e = e + 2^ e^ 4- B, and A = — ^ — B ; 

whence 

^ = G j 1 4- sin2 I {^— m - e + jB ~ ~ m ej — sin^ / cos^ ly-^me—-^ + 3 AjL 
which is the same expression as that obtained by Mr. Airy in the Philosophical Trans- 
actions, 1826, except that instead of e or , the symbol is used to represent -y- . 

24. The circumstance of the terms arising from the differentiation of V with respect 
to a vanishing, affords an easy method of extending Clairaut's theorem indefi- 
nitely, without calculating the value of V. 

It may be shown independently, that these terms cancel each other in all cases. 
The («4- l)th term of V for the portion including the point is 

The corresponding tei^m for the other portion is 
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and it is evident, by inspection of these functions, that that portion of the sum of their 
partial differential coefficients, which arises from differentiating with respect to r 
under the signs of integration, is equal to zero. This being the case, it is not neces- 
sary to know the forms of the functions M a, N a, P a, &c., nor their numerical coeffi- 
cients; but only the functions of ^, by which they are respectively multiplied ; and 
these are Laplace's coefficients. 

Thus the equation of equilibrium at the surface would be 



1 / „ l\Na , 1 / , 6 , . 3\Pa / « 15 ^ , 5 ^, 5 \ Q« 



+ g- (1 — /^2) r2 m G, where - =~(1 -fs|«'2 + A/K,4-f-D ^^), suppose. 

Expand r, recollecting that N a, P a, and Q a are of the 1st, 2nd, and 3rd orders 
respectively, and we have three equations to determine g, A and D. 

By differentiating C with respect to r, and eliminating N a, P «, and Q a by these 
three equations, we have the resolved force in r, which divided by the cosine of the 
angle of the vertical gives g exactly as in (23.). 

It is evident that this may be carried on indefinitely ; and to any order, without 
finding g for the next lower order. 
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